High-resolution Digital Surface Model of the 2021 eruption deposit of Cumbre Vieja volcano, La Palma, Spain

Identifying accurate topographic variations associated with volcanic eruptions plays a key role in obtaining information on eruptive parameters, volcano structure, input data for volcano processes modelling, and civil protection and recovery actions. The 2021 eruption of Cumbre Vieja volcano is the largest eruptive event in the recorded history for La Palma Island. Over the course of almost 3 months, the volcano produced profound morphological changes in the landscape affecting both the natural and the anthropic environment over an area of tens of km2. We present the results of a UAS (Unoccupied Aircraft System) survey consisting of >12,000 photographs coupled with Structure-from-Motion photogrammetry that allowed us to produce a very-high-resolution (0.2 m/pixel) Digital Surface Model (DSM). We characterised the surface topography of the newly formed volcanic landforms and produced an elevation difference map by differencing our survey and a pre-event surface, identifying morphological changes in detail. The present DSM, the first one with such a high resolution to our knowledge, represents a relevant contribution to both the scientific community and the local authorities.


Background & Summary
The morphology of active volcanoes is dynamically shaped by eruptive activity and erosional processes acting at different timescales. Consequently, a precise digital elevation model is fundamental for mapping volcanic hazards, modelling volcanic processes, and complementing further analysis. Furthermore, in urbanised areas, detailed post-eruption topography is important for land recovery actions. Volcano morphologies can be quantified using different techniques [1][2][3][4][5][6][7][8][9] . Recently, the increased capability of UASs and their applications for aerial observation 10,11 , together with the parallel development of Structure-from-Motion (SfM) process 12 , brought important and valuable advantages compared to the classical ground-based, satellite, and crewed aircraft surveys. Nowadays, UAS-based photogrammetry is routinely applied on volcanoes to obtain very-high-resolution DSMs [13][14][15][16] .
Cumbre Vieja is the active volcanic rift on La Palma and has seen the largest number of eruptions of the Canary archipelago in historic times 17 , and its 2021 eruption was the largest eruptive event in recorded history for La Palma. The previous eruption occurred in the southern part of the island between September and November 1971. The 2021 eruption was preceded by an unrest phase characterised by increased ground deformation starting from 2009 18 , increased seismicity from 2017 [19][20][21][22] , and detection of geochemical anomalies from 2010 23 . A dramatic evolution of the seismicity began on 11 September 2021 with a seismic swarm characterised by an upward migration of the hypocentres reflecting the rising of magma towards the surface.
The volcanic eruption at Cumbre Vieja started on September 19, interrupting its 50-years-long period of quiescence, and lasted until December 13 (85 days and 8 hours 24 ). During this period, the volcanic activity was distributed along a fissure where a multiple-vents volcanic edifice formed (called "Volcán de Tajogaite''). The explosive activity was characterised by alternating strombolian explosions and lava fountaining episodes, accompanied by abundant lava effusion. All such phenomena produced profound morphological changes in the landscape and severely affected settlements and industry. A total of about 12 km 2 of territory, more than 1,600 buildings and 200 hectares of banana plantations (the island's main economic resource after tourism), and important infrastructures (roads, powerlines, waterlines, etc.) were buried and destroyed by lava flows in their 6-km-long-path to the ocean. Here they expanded into two lava deltas, forming new land. In addition, the tephra fallout further affected the whole island and, to a smaller extent, the nearby islands of El Hierro, La Gomera and Tenerife.
Here, we present the results of a UAS survey carried out between 24 and 28 January 2022 using a DJI Phantom 4 RTK (real-time kinematic). The aerial images were georeferenced using an onboard RTK receiver capable of cm-level positioning accuracy. The dataset was then processed using Structure-from-Motion (SfM) photogrammetry and 40 Ground Control Points (GCPs) acquired between 23 and 27 January using the Differential Global Navigation Satellite System (DGNSS) positioning. This allowed us to achieve horizontal and vertical centimetre accuracy and to produce a very-high-resolution (0.2 m/pixel) Digital Surface Model (DSM) and orthophotomosaic (0.1 m/pixel), covering an area of about 17 km 2 .
Topographic change detection was obtained by differencing our survey and a pre-event (2015) 2m-pixel DTM 25 , thus identifying elevation changes at decimetre level precision. We characterised the whole topography of the new volcanic edifice and related lava field to detect elevation, areal, and volume variations.
Summary of the main findings: • Subaerial deposit of lava flows and proximal fallout: volume 217.4 ± 6.6 Mm 3 (voids in the lava field and submerged portion of the two deltas are not considered), subaerial deposition rate 29.5 m 3 /s. In a previous survey carried out on 27/9, 35.8 ± 3.0 Mm 3 and 59.2 m 3 /s were recorded 26 . Moreover, considering the volume difference between the 27/09 survey 26  The present DSM represents a relevant contribution to both the scientific community and the authorities in charge of the restoration activities management.

UAS survey and DSM generation.
We conducted a photographic survey campaign ( Fig. 1 and Table 1) between 24 and 28 January 2022, collecting multiple sets of UAS-based high-resolution imagery. We acquired over 12000 aerial pictures using a DJI Phantom 4 RTK UAS with a 1" CMOS 20MP and a field of view (FOV) 84° and 8,8 mm/24 mm (35 mm equivalent) focal length lens.
A total of 10 multi-flight missions were conducted for the survey, for a cumulative flight path of over 800 km (Fig. 1). All flights except two were nadir image data collection missions, conducted at an approximate altitude of 200 metres above ground level (a.g.l.), resulting in a nominal ground-sampling-distance (GSD) of 5.4 centimetres per-pixel. The 24 and 28 January 2022 flights carried out in the area of the cone were both nadir and oblique image data collection missions conducted at a variable altitude of 50-200 metres a.g.l. For the nadir flights, we flew the UAS using predefined missions. Flight planning was designed with 80% forward and side overlap at ground level. Before each flight, we adjusted the camera's digital ISO, aperture, and shutter speed according to ambient light conditions.
With respect to other terrains, several additional difficulties characterised the aerial photographic survey campaign at Cumbre Vieja. The cone area has a highly irregular topography, characterised by notched craters and slopes. In addition, viewing conditions at Cumbre Vieja were still partially limited by the presence of vapour/gas plumes and, at times, by atmospheric haze and clouds.
The data on camera position were collected using GNSS-RTK information embedded in the image metadata by means of a DJI D-RTK 2 Mobile Station. In addition, 40 ground control points (GCPs) were distributed along the outer boundary of the lava flow and in the cone area to establish survey control (Fig. 1). In detail, 33 points were used as proper GCPs (i.e., used to georeference and scale the photogrammetric model and for camera calibration purposes), whereas 7 points were used as checkpoints (i.e., not directly used in the photogrammetric modelling process but available to check the accuracy of the generated model). GCPs were measured with a GNSS survey using a DJI D-RTK 2 Mobile Station in real-time kinematic (RTK) mode, with differential corrections sent in real-time by the Instituto Geográfico Nacional differential positioning service available at https:// www.ign.es/web/ign/portal/gds-gnss-tiempo-real. The surveyed GCPs have an accuracy of 1-2 cm in horizontal coordinates and 2-4 cm in elevation.
Following image collection, we culled the photoset, removing dark and/or blurry photos. We then pro- www.nature.com/scientificdata www.nature.com/scientificdata/ cloud by using the "filter by confidence" tool and by manually removing anomalous floating points caused by the presence of the volcanic plume; and (6) generation of DSMs and orthomosaics. We set the processing parameters in Agisoft Metashape ® to "high" for photo alignment accuracy and "high" quality and "aggressive" depth filtering for dense point cloud generation. For the details of the photogrammetric survey data and elaboration refer to Table 1.
We generated a 0.2 m/pixel DSM (Fig. 2) and a 0.1 m/pixel orthophotomosaic, covering an area of about 17 km 2 . Unlike a digital elevation model (DEM), the DSM represents the elevation of the highest object within the bounds of a cell. Vegetation, buildings, and other objects have not been removed from the data.
The dataset was processed and analysed in the REGCAN95/UTM zone 28 N [EPSG:4083] Coordinate System. The transformation from ellipsoidal to orthometric heights has been performed using the Geoid model EGM08-REDNAP (https://datos-geodesia.ign.es/geoide/). (Fig. 3) was obtained by differencing our surveys and a pre-event 2m-pixel DTM acquired in 2015 for the Spanish PNOA-LiDAR project 25 . The assumption is that between the acquisition of the pre-event DTM (2015) and the beginning of the volcanic activity (19 September 2021), no significant height variation took place in the study area, so that elevation differences obtained in our analysis are mainly linked to the volcanic eruption. To subtract the post-and pre-event surveys, we resampled our DSM to 2 m/pixel resolution (same resolution as the 2015 DTM). Considering the vertical Root Mean Square  www.nature.com/scientificdata www.nature.com/scientificdata/ Error (RMSE) of 0.26 m for our model (before resampling), we set the threshold elevation change (minimum level of detection or minimum elevation change that can confidently be considered a true change) to 0.5 m. It is worth mentioning that the pre-event reference surface is a DTM while our product is a DSM. Such difference must be considered when subtracting both layers as height contributions from vegetation and buildings are still present on the DSM. However, the contribution of such areas is negligible as they are not present above the lava flows and in the cone area.

Data Records
The data record consists of a high-resolution (0.2 m/pixel) photogrammetric Digital Surface Model processed from survey campaign photographs using Agisoft Metashape ® . Details of the photogrammetric survey data and elaboration are summarised in Table 1

technical Validation
Errors in our photogrammetrically-generated DSM result from a complex interplay of geometric and physical parameters, such as image scale, GSD, camera network geometry (nadiral, cross, oblique strips), percentages of image overlap (forward and sidelap), camera shutter speed and exposure settings, lens specifications, image sharpness, camera calibrations, flight design (e.g., flight-line geometry and altitude), surface texture and albedo, lighting conditions, accuracy and distribution of GCPs, disturbances from volcanic activity, as well as on processing: SfM, BBA, image matching, point cloud noise, and outlier removal algorithms.
We therefore applied several strategies to mitigate errors, among which the most important were the following: (1) the use of fast (>1/400 s) camera shutter speeds (i.e., exposure times) whenever possible, (2) the variation of flight altitudes and camera orientation, (3) the application of best practices for processing in Agisoft Metashape, (e.g. 12 ), and (4) the removal of sparse cloud points with large uncertainty via Metashape's gradual selection tools.
The technical quality of the reconstructed DSM was assessed by using the survey report generated by Agisoft Metashape ® and by comparing our DSM to a pre-event DTM 25 . According to the Agisoft Metashape ® survey report the GCPs and check points error estimates are as follows: the total GCPs Root Mean Square Error (RMSE) www.nature.com/scientificdata www.nature.com/scientificdata/ is 6.18 cm and the total check points RMSE is 14.59 cm ( Table 2). The residual elevation difference with respect to a 2 m/pixel pre-event (2015) DTM 25 extracted at 13 check points placed in the unchanged regions of our DSM was used as an additional indication of the vertical RMSE, which is 0.26 m. Our model is thus sufficiently accurate for the scale of changes reported in this study.

Code availability
No custom code was used to generate or process the data described in the manuscript.   Table 2. Ground control points (GCPs) and check points Root Mean Square Error (RMSE). X -Easting, Y -Northing, Z -Altitude.